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Abstract. A bond-disordered two-dimensional Ising model is used to simulate Kauzmann's mechanism 
of vitrification in liquids, by a Glauber Monte Carlo simulation. The rearrangement of configurations is 
achieved by allowing impurity bonds to hop to nearest neighbors at the same rate as the spins flip. For 
slow cooling, the theoretical minimum energy configuration is approached, characterized by an amorphous 
distribution of locally optimally arranged impurity bonds. Rapid cooling to low temperatures regularly 
finds bond configurations of higher energy, which are both a priori rare and severely restrictive to spin 
movement, providing a simple realization of kinetic vitrification. A supercooled liquid regime is also found, 
and characterized by a change in sign of the field derivative of the spin-glass susceptibility at a finite 
temperature. 

PACS. 64.70.Pf Glass transitions - 75.10.Nr Spin-glass and other random models 



1 Introduction 

In a seminal article maiiy years ago, Kauzmann pro- 
vided the now standard [g physical picture of vitrification 
in liquids, which emphasized that obstruction of config- 
urational rearrangement by molecular interactions could 
lead to such an inefficient sampling of phase space, that 
thermodynamic equilibrium could not be reached for very 
long times. The ability to change configuration practically 
disappears at the glass transition, while orientational dis- 
order does not change significantly, being high and rela- 
tively static in both the supercooled liquid and the glass. 
For our simple purposes, the scenario may be compared 
to a log-jam. There is nothing in the Hamiltonian to pre- 
vent orderly stacking, and indeed, vitreous liquids possess 
a crystal state of lower free energy. But once the logs are 
jammed, it is (very) improbable they will be ordered by 
random pushing and pulling. 

By contrast, in spin glasses jsj^, the Hamiltonian it- 
self forbids crystal order, because of impurities, or lat- 
tice topology. Now disorder in terms of the original mi- 
croscopic variables becomes an equilibrium property at 
all temperatures, leaving open the possibility of a true 
thermodynamic phase transition in terms of some other, 
less-than-obvious order parameter. In searching for it, the 
basic concept is frustration, introduced ||^ to distinguish 
'intrinsically' disordered states from those where apparent 
disorder can in fact be removed by redefining parameters 
of the Hamiltonian. These 'serious' and 'trivial' kinds of 
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disorder may be called 'glassy' and 'amorphous', respec- 
tively. 

In the present article, we bring the two subjects to- 
gether, to obtain a toy model of kinetic vitrification. In 
brief, our algorithm gives a structural glass in a spin- 
glass setting. We begin with the discrete bond-disordered 
Ising model in two dimensions. It is known [|| that it 
has no spin-glass transition at T > 0. Topologically frus- 
trated configurations are too rare in two dimensions to 
contribute significantly to the sample-averaged free en- 
ergy, when each random arrangement of impurities con- 
tributes a sample with equal a priori probability. We refer 
to this standard quenched-disorder picture as 'frozen' or 
'frozen bonds', reserving the word 'quenched' as follows. 

The model we consider is obtained from the above by 
allowing the impurity bonds to move. Throughout the ar- 
ticle, 'quenching' means that this system, with movable 
bonds, is first equilibrated at high temperature, and the 
temperature then brought down to some low value in a 
single step, after which the evolution continues at that 
low temperature. Since the bonds move, there is no aver- 
aging over samples; the complete impurity configuration 
space is available to a single sample. Of course, the effi- 
ciency with which this space is explored after the quench 
critically depends on the nature of bond movement. We 
focus on the case where they move locally, i. e. a bond can 
only exchange places with another bond, impinging on the 
same site. Our main point is that such an algorithm, which 
would be very inconvenient to study the ground state of 
the disordered Ising model, because it is easily trapped in 
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metastable states, is very convenient from the structural- 
glass point of view, because the metastable states in which 
it does get trapped are quite exceptional among all possi- 
ble impurity arrangements, and are analogous to vitreous 
states in the structural sense. 

It is easy to imagine how this happens: the bonds 
keep moving until they get stuck. These 'stuck' config- 
urations are precisely the vitreous ones, which are (by 
definition) those in which evolution appears to cease be- 
fore the ground state is reached. Various signals of their 
exceptional nature are studied throughout the article; let 
us only mention that their energy is markedly lower than 
the average energy in the frozen scenario, in fact close 
to the theoretical ground state energy (which can also be 
reached, by slow cooling). Nevertheless, it is finitely above 
the ground state, making them metastable. 

In the present work, we consider only the fast evolution 
after a quench, leading to the above-mentioned metastable 
states, and the evolution of these states at intermediate 
times, but not their relaxation to the ground state, occur- 
ing at much longer time scales. The ground state (annealed 
system) is a comparatively trivial amorphous configura- 
tion, and our only interest in it is that it establishes the 
formal metastability of the long-lived quenched states. In 
this, of course, our toy model is different from a true vit- 
reous liquid, where the ground state is an ordered crystal. 
The model analogue of melting involves no order-disorder 
transition. 

In the description of our results, we use established 
language {e.g., melting) and theoretical tools {e.g., the 
spin-glass susceptibility) somewhat outside their original 
context. The model is purely kinetic, without a thermo- 
dynamic phase transition, so the interesting (vitreous) be- 
haviour appears in quasi-equilibrium, established at the 
above-mentioned intermediate time scale. We are guided 
by an intuitive analogy between spins and bonds in the 
disordered Ising model, on the one side, and the orienta- 
tional and configurational degrees of freedom, respectively, 
of a glass-former, on the other. 

One class of Ising-like models, often used to study vit- 
rification, are constrained kinetics models, such as the 
Frederickson- Andersen model 0. They have no intrinsic 
disorder, but frustration is introduced by dynamical re- 
strictions on changing the spin states, which take the part 
of configurational variables. Among work on variants of 
the disordered Ising model, our program is closest to that 
of the frustrated percolation problem ||,|^. However, that 
approach views impurity bonds as obstacles to rearrange- 
ment, playing a role similar to the dynamical restrictions, 
mentioned above. In our work, bonds are imagined to rep- 
resent the configurational degree of freedom itself, so their 
movement becomes of physical interest. It is responsible 
for the appearance of metastable states, the 'supercooled 
bond liquid', which have a distinct physical characteris- 
tic, a first-order zero in the field derivative of the spin 
glass susceptibility, not shared by the states on the cool- 
ing curve connected to the ground-state manifold. This 
liquid, and the accompanying glass, are the main subject 
of the present article. 



2 Model and algorithm 

The model is the bond-disordered two-dimensional Ising 
model on a square lattice p^ : 

H=-Y. J^JS^S, - BJ2S^^ -J2 -B-M, 

iv) i (ij) 

(1) 

where the 'spins' take the values Si = ±1, as well the 
bonds, Jij = ±1. Here M is the total magnetization, cal- 
culated as above, not an independent parameter. The frac- 
tion of antiferromagnetic (negative Jij) bonds is a fixed 
number p; in this paper, p = 70%. The key new feature 
is the ability of bonds to move; this is considered to be 
purely diffusive, so there is no kinetic contribution to the 
Hamiltonian. 

Spin updates proceed by flipping a spin at random; 
the criterion for acceptance is Glauber's. Between spin 
updates are bond updates; the test move is to exchange 
two bonds of the opposite sign, impinging on the same 
site. This is accepted or rejected, again by the Glauber 
criterion. Which bonds are exchanged, and on which site, 
is chosen at random, taking care, of course, that the two 
bonds involved are indeed of the opposite sign. There is 
no regular 'sweeping' of the lattice, to avoid artificial cor- 
relations between spin and bond updates. 

There still remains a physical freedom in the update 
protocol. One might consider, in general, flipping s spins 
and moving b bonds, before submitting the new configura- 
tion for acceptance. These are large steps in the random- 
walker picture. They are unlikely to be accepted in gen- 
eral, but once the system is trapped in a glassy state, such 
steps may significantly affect its chance to get out. These 
rare 'cooperative' fluctuations are considered responsible 
for re-crystallization of vitreous liquids in the Kauzmann 
picture, sometimes after 'archeologically' long times, as in 
window glass. 

In this work, we only focus on driving the system into 
a glassy state, not getting it out. On this 'laboratory' time 
scale, alternating one-spin and one-bond updates, each 
subject to acceptance in turn, provide the dominant quasi- 
equilibration mechanism. The very-long-time behavior of 
the model is the subject of a limited investigation in Sec- 
tion ^ below, but still neglecting multistep (cooperative) 
processes, which could be important for the asymptotic 
regime. Parenthetically, we note that mode-coupling the- 
ory introduces an artificial thermodynamic transition, and 
therefore certainly has the wrong asymptotics, but even 
this does not affect the intermediate time scales, for which 
it was devised . 

We consider two situations: annealed, where the sys- 
tem is equilibrated at high temperature (T = 5, but any 
T > 2 will do), and then slowly cooled (in steps of 0.05, 
waiting to equilibrate between temperature steps); and 
quenched, where it is brought down to the desired temper- 
ature in a single step from the high-temperature equilib- 
rium, and then quasi-equilibrated. The high-temperature 
starting point is generated anew for each quenched point, 
so the agreement of annealed and quenched curves for not- 
too-low temperatures is also a check of equilibration. Sim- 
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ilarly, the fact, that the annealed curve nearly reaches the 
ideal (defect-free) ground state energy at T = 0.05, our 
lowest temperature, is indication that the cooling rate for 
annealment {AT 10~^ per MC step) is sufficiently slow 
for our purposes. 

For comparison, we also simulate some systems with 
frozen bonds. This opens the question of sample-averaging. 
We find that self-averaging is sufficient for our limited pur- 
pose, since the samples involved are rather large (100 x 100 
and larger). To be safe, we have made some runs, averaging 
batches of up to 500 small (~ 10 x 10) samples, and in- 
deed no claims made in the paper, based on the large self- 
averaged ones, were affected. Note, also, that the frozen 
bond arrangement at high temperature is different for each 
cooling run, so the above-mentioned check of equilibration 
is also a check of self-averaging in the frozen case. 



3 Main result 




Fig. 1. Energy per site for p = 0.7 and B — 1. Full curves: 
annealed. Dashed curves; suddenly cooled. Lower (upper) two 
curves: with (without) bond movement. Open circles; quench 
to T = 0.4 and slow cooling. Symbols at T = 0.05: see Figs. ^ 
and|. C/min = -2.3 [Eq. (|p], lattice 100 x 100. 



All simulations were performed in a uniform external field, 
whose non-trivial range is < B < 4. Then the theoretical 
minimum energy per site of the annealed system with a 
fraction p of antiferromagnetic bonds is 



C/mi„ = -2-(l-p)B, 



(2) 



which assumes the absence of defects: each spin, turned 
opposite to the field, is connected to its neighboring sites 
by antiferromagnetic bonds, which are all used up in this 
way. These 'crosses' carry no frustration, obviously: a si- 
multaneous sign change of the spin and its four bonds 
leads to a perfect ferromagnet. There is a ground-state 
entropy, associated with their random placing on the lat- 
tice, so this ideal state is amorphous. It is our model's 
analogue of the crystal state of vitreous liquids. 

A convenient robust measure of additional disorder 
is the energy mismatch between this reference state and 
one obtained in an actual run. This can be extended to 
any temperature by using the reference system's cooling 
curves, once it is checked that its zero-temperature limit 
is near the theoretical amorphous one, i.e., that it is suf- 
ficiently annealed. 

Figure ^ shows the energies per site for a single cool- 
ing run, in an external field B = 1 and p = 0.7. Four 
curves are given, annealed and suddenly cooled, with and 
without bond movement allowed. Two energy scales ap- 
pear: a higher one, where the curves with bond move- 
ment separate, and a lower one, distinguishing annealed 
and suddenly cooled systems. Unexpectedly, the higher 
scale is about the square root of the coordination, which 
is the mean- field transition temperature of the model (|l|) , 
Tmf = \/^ = 2; this may be a coincidence. The lower one, 
denoted Tq, is the model 's m elting temperature, a notion 
made precise in section 4.1 below. The main interest in 
this work is the range beneath it. 

Below To , the dashed curves rise with falling tempera- 
ture, because each point on them is generated anew by a 
drop from high temperature. The lower the temperature 
to which they are quenched, the 'worse' (energetically) 



Fig. 2. A 50 X 50 segment of the quenched state, marked 
by a filled square in Fig. 0. Only uncompensated bonds are 
shown. The preponderance of singletons indicates nearly all 
the disorder is glassy. 



the state in which they end. In the next section, we shall 
discuss a third type of cooling run, where the system is 
cooled slowly after a quench from high temperature. An 
example of this is also given in the figure, showing the 
thermodynamically expected behavior. 

While the two dashed curves look similar, they refer 
to very different states. In Figures § and |[ we draw the 
respective lattices at low temperature, presenting only un- 
compensated bonds, i.e. those which are of the 'wrong' 
sign relative to the spins they connect. The exact ground 
states have as few of these as possible (a nonzero num- 
ber may be topologically unavoidable 0). The quenched 
state obviously has more, but the point is that they nearly 
always appear as well-separated singletons, carrying topo- 
logical frustration. In other words, the frustration of our 
vitreous states is near the maximal possible one, for a 
given energy difference from the ground state. In conse- 
quence, spin orientation is fixed by a 'majority vote' of 
bonds at every single site (almost). The comparatively 
small energy difference between the ground and vitreous 
states is misleading; the two are separated by a topologi- 
cal barrier. The fact that such a priori rare states can be 
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Fig. 3. A 50 X 50 segment of the frozen, rapidly cooled state, 
marked by a filled triangle in Fig. ^. Only uncompensated 
bonds are shown. Much of the disorder is amorphous. 
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generated with so little effort, about 50 updates per site, 
is an important practical point of this work. 

By contrast, in the case with frozen bonds, there ap- 
pear many sites with two compensated and two uncom- 
pensated bonds, making the spin orientation energetically 
indifferent. Many of the corresponding plaquettes are not 
frustrated, so the energy difference relative to the ground 
state is largely due to trivial (amorphous) disorder. Be- 
cause of this, the magnetization of the frozen-bond system 
is a sum of fluctuating contributions even at low temper- 
ature, making it constant in time only on the average, 
while the magnetization of the vitreous (quenched) sys- 
tem is a sum of spins which are individually fixed in time 
for the lifetime of the glass, as would be the case for all 
times with a perfect Ising ferromagnet at sufficiently low 
temperature. 

Figure ^ shows the magnetization. For a frozen bond 
configuration, it barely distinguishes between annealing 
and quenching. What is unexpected is the difference be- 
tween the annealed and quenched magnetizations, when 
bond movement is allowed. The annealed magnetization 
at T w is higher than when bonds are frozen, reason- 
ably enough. The quenched one, however, is much lower, 
indicating a system in which spin orientation is controlled 



by an internal field. It is very sensitive to the Tq scale of 
Fig. 1^, while the same scale is hardly visible in the other 
three curves. 

It remains to show that the state quenched (in our 
sense, with mobile bonds) below Tg « 0.2 fulfills the cri- 
terion for a spin glass. We call it a 'bond glass', to em- 
phasize the role of bond kinetics in its formation; it is 
not in thermodynamic equilibrium. The transition region 
Tg w 0.2 < T < 0.8 w To in Figure | is then a 'supercooled 
bond liquid', characterized by a scale Ti « 0.5. For tem- 
peratures above Tq , we have a bond liquid. CarefuU cooling 
below the melting temperature Tq leads to the amorphous 
solid, described above. (Note the small downturn below Tq 
in the slowly cooled curve, which seemed headed towards 
the theoretical maximum magnetization M — 0.3. It is due 
to defects in the amorphous structure. The unavoidability 
of these is the subject of bond-percolation studies [|Tl|].) 

We consider the spin-glass susceptibility, which is a 
time-asymptotic four-spin correlation function 



XSG = V q^P{q)dq, 



(3) 



where the 'volume' V is the number of spins in the lattice, 
and P{q) is the distribution of the quantity 



Q ~ lim — — 



1 
At 



At 

E 



i—1 



t)S, 



(2), 



(4) 



where the superscripts on the spins denote two copies of 
the system, independently evolved, and 'infinity' is some 
intermediate time scale, before ergodicity is re-established, 
as discussed below. For each temperature, we first quasi- 
equilibrate the system, as to produce a point in the first 
two figures, then make two copies of the lattice, identi- 
cal in both spins and bonds, so an artificial divergence 
is introduced in XSG at t — t — 0. This is allowed to 
decay by evolving the copies independently for a further 
r = 5 X 10^ MC steps, each consisting of a spin and a bond 
trial, or only a spin trial, for the cases with frozen bonds. 
The distribution P{q) is obtained by considering the next 
At = 2.2 X 10^ steps. 

In the frozen case, this is done differently : one ran- 
domizes the spins in both copies of the lattice, with the 
same bond arrangement, and waits for the correlation 
function to build itself up again from zero. We cannot 
do that, because the bonds also move: once the spins are 
randomized, bonds start flowing, and the lattice copies 
end up as distinct samples, as soon as the vitreous state 
is re-established, which takes a very short '/3-relaxation' 
time, only ~ 50 updates per site in our algorithm. In the 
usual picture of a glass free energy landscape, the two 
copies end in distinct deep and narrow minima, and the 
scale associated with tunneling between such minima is 
the long 'of-relaxation' time. In the same language, the 
time window At is opened after 'inflnity' ~ 20 times the 
/3-relaxation period, i.e. in the quasi-equilibrated regime, 
Tfj < T < Ta for all T, which is the regime of interest. 

The use of the spin-glass susceptibility should not be 
construed to mean, that there is a true thermodynamic 
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Fig. 5. Normalized spin-glass susceptibility on a 120 x 120 
lattice, for a number of quenching scenarios. Empty symbols: 
B = 1. Filled symbols: B = 0'^ . Circles: quench with movable 
bonds. Squares: same as circles, but no bond movement after 
r = t = 0. Thick (thin) dashed line: B = 0+ (S = 1), random 
frozen bonds. Zero-field Tg, Ti and To are also given. 

transition in the model. As long as the bonds move, all the 
reported phenomena are kinetic in origin. Asymptotically, 
XSG is zero at all temperatures, but the time needed to 
see this increases exponentially with falling temperature, 
as discussed in the next section. 

The results are shown in Figure ^, which we discuss for 
the remainder of this section. The open circles correspond 
to the quenched bonds in Figure ^. The range T < 0.2, 
where the magnetization saturates at M « 0.15, is clearly 
linked to a divergence in xsG- The crossover region, 0.2 < 
T < 0.8, is also evident in both figures. 

The thin dashed line corresponds to the frozen-bond 
scenario, whose magnetization is given by the nearly hori- 
zontal dashed line in Figure ^. It saturates at Xsg/V = 1 
as well, which may cause doubt about the alleged special 
nature of the state, found by bond movement. To dis- 
pel it, we repeat all of the above with the magnetic field 
turned off. (We put B = 0.001 to break inversion symme- 
try, and denote this B = 0+.) The frozen-bond suscepti- 
bility then drops to the lower, thick dashed curve, which 
no longer shows saturation. This indicates its divergence 
at finite field was in fact trivial, due to the difficulty of 
fiipping spins when T <C -B, so the spins on the two copies 
'marched together' for a long time. 

The quenched susceptibility shows precisely the oppo- 
site behavior, given by the filled circles. For B = 0+ it 
rises above the 5 = 1 points, to achieve nearly perfect 
divergence for T < 0.3, but drops below them in the high- 
temperature part, so the net effect is a narrower crossover 
(supercooled liquid) region, Tg = 0.25 < T < 0.75 = To, 
its center (half-point of xso/y) remaining at = 0.5. 
The low-temperature divergence in the quenched case is 
due to the bonds 'holding' the spins, in place of any exter- 
nal field. The field widens the supercooled liquid regime, 
and produces a less perfect glass (slightly smaller xsG be- 
low T = 0.2). At B = nothing is visible in the magne- 
tization, which fluctuates around zero at all temperatures. 

We now turn to the role of bond rearrangement in the 
supercooled liquid. We restart the evolution from the same 



C(0,T) 




10° io' id'' 10^ 10* lo'^ 10° io' 



T [MCS/site] 

Fig. 6. Spin autocorrelation function, quenched case, high- 
temperature initial state. Left to right by arrows: T = 0.4 
(L = 25), T = 0.3 (L = 15), T = 0.25 (L = 15). Arrows: 
dip marking Td- White squares: points with abscissa given by 
Eq. (^. Dashed line: a minimal r^, for all cases. 

copies of the lattice as gave the circles in Figure ^, but now 
disallow further bond movement. In other words, we pro- 
ceed as with frozen bonds (the dashed curves) , but for the 
special bond configurations found by quenching (Fig. 
not random ones (Fig. ^). Then the open and filled circles 
go over into the respective squares. The supercooled liquid 
has disappeared, replaced by a glass up to T « T^. This 
is in accord with the idea, that changing configurations is 
essential for the liquid. It may appear that the crossover 
has merely moved to higher temperatures, not shown in 
the figure. But T = 0.8 is already the melting tempera- 
ture in Figure ^, so the system cannot be 'supercooled' at 
T > Tq. It melts, as it eventually must, by a trivial high- 
temperature degradation of the initial spin orientations. 
However, the rarity of the bond configurations involved 
is evident in the lack of field-dependence of the suscepti- 
bility (the squares nearly overlap) , in sharp contrast with 
the behavior of the dashed curves. 



4 Equilibration 

The basic point of this article is that the algorithm with 
bond movement is physically interesting in itself, as a 
model realization of a structural glass. At present, we have 
no physical argument that the degrees of freedom of any 
particular real glass map onto our model. So one has at 
best an analogy, or 'toy model', and the idea is to use 
the analogy to guide model investigations. Some quanti- 
tative evidence will be presented in this section, that the 
properties uncovered in this way are not spurious. 



4.1 Short and long times 

In Figure ^ we show the spin autocorrelation function 

cit,T)^-Y,{s,mit+T)), (5) 
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where the brackets denote an average over time-windows 
of 10 MCS/site. The reference state {t — 0) is taken to 
be the high-temperature equihbrium (T — 5) before the 
quench. A typical two-time evolution is observed. After 
a short /3-relaxation time rp, 50-100 MCS/site, the sys- 
tem enters a long-lived low-temperature state. This inter- 
mediate period, in which evolution almost seems to stop, 
ends at a time Tqi, when C(0,t) begins to decay further, 
reaching zero only after a much longer time Tq,2. The lat- 
ter is usually called Tq, the scale at which ergodicity is 
re-established. This is the time needed to tunnel between 
the 'deep and narrow' minima characterizing the glass free 
energy landscape. In the same picture, Tai would be the 
time needed to exit the original minimum, which trapped 
the system at t k, rp. 

It is easy to estimate Ta2 in the present model, be- 
cause Fig. |2| indicates maximal topological frustration is 
a reasonable hypothesis. This implies, roughly, that each 
spin needs to flip against the local field before the state is 
completely dissolved. In other words, the reciprocal prob- 
ability for a single spin to flip is the waiting time for disso- 
lution, measured in updates per site. The energy difference 
for a flip is AE = 4 or AE = 8, depending on whether the 
spin is held in place by three or four compensated bonds. 
Taking the lesser of these gives an underestimate. 



(6) 



C(t,T) 



exp ( — ) < Ta2 [MCS/site] 



It 



which we find to be logarithmically correct (Fig. [ 
predicts Ta2 > 5 x lO'^'* updates per site at T = 0.05. 

The success of the simple estimate ^ indicates that 
the algorithm with nearest-neighbor bond diffusion real- 
izes, more or less, the maximal time scales possible in the 
model, for a purely kinetic scenario. But it also implies 
that the scenario is indeed kinetic: nothing really happens 
at Tg, except that a time scale becomes too long to follow. 
This is the reason the present approach is computationally 
so cheap. One knows in advance that, say, XSG in Fig. o 
is asymptotically zero. The question is, for what delay [r 
in Eq. (^] it is most informative of the processes in the 
system. This turns out to be fairly modest, ~ 10^ updates 
per site. By then, the value of Tg is as well established as 
for 10^ updates, while the Tq for most of the higher tem- 
peratures are still not reached, giving the largest possible 
range in temperature, for which < XSG < 1. Of course, 
there is a highest temperature, where Ta2 — > t/s ~ 50. This 
is just the melting temperature T^: exp(4/To) ~ 50 gives 
To « 1, as found in Fig. |l|. Above it, xsG is zero at all 
times. 



4.2 Aging 

After the initial relaxation has passed, the supercooled 
liquid/glass spends a long time in a state which is nomi- 
nally metastable, since its energy is above that of the an- 
nealed state (Fig. However, it is equilibrated from the 
point of view of correlations; there is no 'aging' [|l^ in the 
metastable state. In such a quasi-equilibrium, the autocor- 
relation function (H) does not depend on choice of origin t. 




t [MCS/site] 

Fig. 7. Spin autocorrelation function, quenched case, for t > 
Ti3. Bottom to top: 1) T = 0.6, i = 25, r = 10^ MCS/site; 
2)T = 0.6, L = 15, r = 10^ MCS/site; 3) T = 0.2, both cases. 
Time-averaging is over 10* updates per site. 



apart from noise, as shown in Figure 0. Here all t > r/3, so 
the initial (reference) state is always a quasi-equilibrated 
one. The amplitude of the noise generally increases with r 
and the temperature, while it decreases as the system gets 
larger. At T = 0.2, C(<, r) = 1, and there is no noise, for 
the trivial reason, that the state below Tg does not evolve 
at the scale of the waiting times exhibited here. 



4.3 Rattling and relaxation 

It is illustrative of the qualitative behavior of the model, 
to look at the acceptance rates of various kinds of trial 
moves. They are the inverses of the time scales involved. 
We introduce a third type of cooling run, most similar to a 
real experiment. In zero field, the system is first quenched 
into the supercooled region (below Tq), and then slowly 
cooled through Tg, in steps of 0.02. We look at the time 
evolution leading to quasi-equilibration (not after it, as in 
the study oixsc)- After the temperature is brought to the 
desired value, a 120 x 120 lattice is subject to 7.2 x 10^ spin 
trials, interlaced with the same number of bond trials. The 
principal quantity monitored is the number of accepted 
relaxation steps, i.e. those in which the energy changes 
either way (positive changes are counted as relaxation, 
because they also mean the system is 'working'). 

The estimate (^) indicates that by lowering the tem- 
perature below Tg , one can really prevent the system from 
reaching the annealed state in any fixed time. In addi- 
tion to the disappearance of spin flips, the nature of bond 
movement changes below Tg, and only 'rattling' |^ of the 
kind depicted in Figure ^ survives. For our purposes, we 
do not define 'rattling' geometrically, but by the property 
that the move is energetically indifferent, so that it cannot 
contribute to annealing, by definition. 

In Figure |9[ we show the total number of relaxation 
steps, accepted during the run above. Bond relaxation dis- 
appears around Tg, while spin relaxation stops earlier, and 
is a rough order of magnitude less significant everywhere. 
Spin rattling vanishes along with bond relaxation. Thus 
both the energy and magnetization cease to evolve below 
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spite the latter being by definition metastable, as it con- 
sists of liigher-energy configurations. This is further illus- 
tration of the long-lived nature of the metastable quasi- 
equilibrium. In the absence of frustration, the difference 
in energy would translate into a larger number of accepted 
steps in the metastable case. This is actually observed at 
T ~ 0.6, but only because the system is quenched from 
T — 5 there. Even when shedding such a large energy, the 
excess steps are almost exclusively due to spins, indicat- 
ing that the bond ('configurational') degree of freedom is 
most affected by frustration. 



Fig. 8. A 'caged' AF bond, free to rotate around the central 
spin. The other bonds (not shown) are all of the correct sign, 
corresponding to the spins they connect. This bond could be 
freed by flipping any of the four spins in the centers of the 
square sides, except the right-hand one. 
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Fig. 9. Number of accepted steps during quasi-equilibration 
of a large system. Symbols: slow cooling after quench to 
T = 0.6; circles: bond moves; squares: spin flips; filled: re- 
laxation; empty: rattling. Full lines: same as symbols, but an- 
nealed. Dashed lines: spin moves with frozen bonds, upper: 
rattling, lower: relaxation. Vertical dotted line: Tg from Fig. |^. 

Tg. This is also visible in Fig. |l|, the open circles there be- 
coming strictly horizontal below Tg. Bond rattling, on the 
other hand, remains high and relatively constant through- 
out, in keeping with the analogy to 'caged' configurational 
movement. 

This is in marked contrast to the behavior of a quen- 
ched, then slowly cooled system with random frozen bonds, 
shown by dashed lines. Spin moves in this case behave 
roughly as bond moves above. Relaxation stops above Tg, 
but spin rattling is remarkably constant down to lowest 
temperatures, as one would expect from the absence of 
a spin glass in 2D, i.e., of an internal field to hold the 
spins in place. On the other hand, spin rattling in the 3D 
Edwards-Anderson model persists even below the spin- 
glasss transition. This is in contrast to the 2D vitreous 
state, in which there are virtually no loose spins at all 
(Fig.|). 

Finally, wc compare the results of a standard annealing 
run, shown by full lines. They are in excellent quantita- 
tive agreement with the supercooled liquid/glass — de- 



5 Discussion 

The main problem with the identification of the glass in 
Section is that XSG is not by itself a measure of frustra- 
tion. It may function as such in a spin-glass setting, where 
one averages over many random bond configurations, so if 
two lattice copies 'march in lockstep' on the average, one 
can safely assume this is not due to the influence of the 
occassional ordered configuration. 

Not so in our case, where the whole idea is that mobile 
bonds find rare configurations. In particular, they find the 
optimal amorphous ones, which show the same persistence 
of the initial divergence in xsGi as would a perfect ferro- 
magnet. The glassy configurations are not far in energy 
(Fig. |l|), and they, too, contain large admixtures of the 
optimal arrangement (Fig. |^) . 

The problem reduces to this: can one find a qualitative 
criterion to distinguish the 'amorphous solid with defects', 
the uppermost of the magnetization curves in Fig. |^, from 
the 'supercooled liquid/glass', the curve which forks down 
from it below Tq7 Luckily, yes. It is provided by the field- 
dependence of xsGj shown in Figure HO for the amorphous 
solid. The filled circles {B — 0~^) show the same behavior 
as in Figure |^, in line with the doubts expressed above. 
But the B — I curve (empty circles) is now entirely above 
the other one; in the annealed case, turning on the field 
increases XSG temperatures. The criterion follows, 

that if the field derivative of the spin-glass susceptibility 

dxsciT) 



is positive at all temperatures, we have a solid. If it is 
negative around zero and changes sign at a finite T = Ti, 
we have a supercooled liquid. 

This criterion makes physical sense. In the amorphous 
solid, the bonds are as good as frozen, since any rattling 
occurs in an optimally arranged background of spins and 
bonds. The spins opposing the field are in a minority (wit- 
ness a finite magnetization), so increasing the field will 
make the system more rigid, increasing xsG at all tem- 
peratures. 

Not so in the supercooled liquid, where the bonds rat- 
tle against a suboptimal background. At low tempera- 
tures, spin flips induced by the field will increase the range 
of bond movement, which in turn has a cascading influence 
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Fig. 10. Same as Figure |5|, but for the annealed state. Note 
the diflterent field dependence. 

(stochastic, of course) on spins in the bond path, degrad- 
ing the correlation between the two copies of the lattice 
(see Fig. H). At higher temperatures, however, there are 
already many spin flips, so this effect is 'played out'; the 
spin majority orchestrated by the field then increases the 
correlation, like in the solid. The temperature at which 
the two tendencies balance is just T^. 

This does not mean that Ti is a thermodynaniical 
(time-independent) temperature, indeed Tg itself is not, 
in any kinetic model. Given that it refers to the liquid 
state, Ti is subject to the same time evolution as XSG in 
that regime, or the correlation function in Fig. ^ Never- 
theless, the above argument implies that the appearance 
of a point, where the zero-field and finite-field suscepti- 
bilities cross, is generic for the supercooled liquid in this 
model, and will always be present for times between 
and Tq.. Since is so short, it can be observed without 
large numerical effort. 

There are two things to note. First, the criterion for a 
supercooled bond liquid implies (mildly) paradoxical be- 
havior: the bonds 'waiting to break out' make the sus- 
ceptibility decrease in response to an applied field, when 
T < Ti. The field has a net disorganizing effect on the 
spins, because it probes the configurational metastability 
of the supercooled liquid. Put simply, the uniform exter- 
nal field is in competition with the stochastic internal field, 
holding the spins in position for T < Tg. 

Second, the glass itself is defined only quantitatively, 
as the manifold of states found on that part of the cooling 
curve of the supercooled liquid, where xsg/V saturates 
at unity. A qualitative characterization of Tg is not to be 
expected, on general grounds. In this context, the switch 
from circles to squares in Figure ^ may be interpreted, 
that if one takes a 'snapshot' of the system below Tg at an 
arbitrary point in time, no analysis of that one snapshot 
can tell, whether it was taken above or below Tg. 

It is easily shown that the maximally frustrated (vit- 
reous) bond arrangements, found by the quenching algo- 
rithm, are a very small subset of all possible ones. For 
a state to be maximally frustrated, uncompensated bonds 
must be singletons, like in Fig. ||, each well separated from 
the others by a 'layer' of compensated ones, effectively 
'using up' q > I bond positions for each singleton. Thus 



uncompensated bonds can choose only among some 2V/q 
positions, less than all 2V available ones. The number of 
distinct states is largest when the number of uncompen- 
sated bonds is one-half the total number of available po- 
sitions, giving at most 

states of this kind. On the other hand, supposing for sim- 
plicity that one has an equal number of positive and neg- 
ative bonds {p = 0.5), the number of ways to place them 
on the lattice without restriction is 

of which the above q-th root is a vanishingly small frac- 
tion, when V is large. Energetically, the vitreous states are 
well separated from the bulk of the arbitrary ones, witness 
Fig. |l|. They correspond to the 'deep and narrow' minima 
of glass theory. In fact, since their energy is determined 
by the number of singletons, Eq. || is also their degener- 
acy, 'exponentially large', as the general picture assumes. 
Since the algorithm is trapped in only one of them during 
any given run (Fig. one has 'breaking of ergodicity' Q 
for times less than Tq,. 

Let us summarize the analogies between our toy model 
and vitreous liquids. Bond and spin movement correspond 
to the configurational and orientational degrees of free- 
dom, respectively. The diffusivity of bond movement re- 
flects the absence of a kinetic energy contribution to the 
glass transition. The supercooled bond liquid turns into 
a glass if configurational rearrangements are artificially 
turned off. Spin disorder changes very little across the 
glass transition, but spin movement ceases there, as well 
as all relaxation. Configurational (bond) rattling persists 
to lowest temperatures. Somewhat above our original am- 
bitions, the zero-field glass transition in the model occurs 
at about one-third the melting temperature, Tg « To/3, 
an old rule of thumb in glass formers j^. 

6 Conclusion 

We have realized a structural spin glass in the two dimen- 
sional bond-disordered Ising model on a square lattice, 
by taking into account only a small subset of impurity ar- 
rangements. This subset was generated by a quenching al- 
gorithm which imitated vitrification in liquids, by allowing 
diffusion of impurity bonds over the lattice. The crossover 
between the glassy and high-temperature (melted) state 
was identified as a supercooled liquid. To distinguish from 
customary usage of the term 'spin glass', and to empha- 
size our glass is kinetic, not thermodynamic, we propose 
the terms 'bond glass' and 'supercooled bond liquid', re- 
spectively. A macroscopic criterion was shown to distin- 
guish the supercooled liquid from the amorphous solid, 
connected to the ground-state manifold by slow cooling, 
which is this model's analogue of the crystal state of vitre- 
ous liquids. The criterion makes it possible to study field 
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and temperature hysteretic effects on the formation of the 
supercooled bond liquid. 
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